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ABSTRACT 


Lithospheric flexure seaward of deep ocean trenches is 
evident in SEASAT altimeter observations of the marine 
geoid. In fact, mechanical models of lithospheric flexure 
can be tested directly on the SEASAT altimeter data. Ve 
have used a simple elastic model for the oceanic lithosphere 
and, after least squares adjustments, have recovered 
estimates of model parameters including Outer Rise (OR) 
amplitude, OR wavelength and effective lithospheric 
thickness . 


Effective lithospheric thickness have been recovered 
for five regions: the Mariana, the Kuril, the Philippine, 
the Aleutian and the Middle America OR. These results 
support the suggestion of Bodlne et al. (1981) that 
effective thickness, T, increases with age of lithosphere in 
approximate accord with the relation T 2 C • a g e ^ ^ where 
C = 4 km*my~l/2. j n f ac t, our altlmetric results agree more 
closely with this relation than do published results based 
on bathymetric data. The close agreement with the 
thickness-age relation suggests that there is no longer any 
need to assume that significant horizontal compression (cf. 
Bodlne et al., 1981) acts across the Kuril, Marianas and 
Izu-Bonin trenches. This thickness-age relation implies 
that flexural strength of the oceanic lithosphere is 
temperature controlled. 
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INTRODUCTION 

The Outer Rise seaward of deep ocean trenches Is 
generally attributed to flexure of the oceanic lithosphere 
prior to subductlon. Various models have been proposed to 
describe this flexure. They all include a mechanically 
strong lithosphere - loaded from the areward side of the 
trench axis ~ overlying a fluid-like asthenosphere . 
Simplest among them is the uniform elastic plate model which 
has been applied to Outer Rises by Hanks (1971), Watts and 
Talwani (1974), Parsons and Molnar (1976), Caldwell et. al. 
(1976) and others. Other models employ a viscous 
(DeBremaecker , 1977 ), a viscoelastic (Melosh, 1978) or an 
elastic - perfectly plastic lithosphere (McAdoo et. al., 
1978). More complex models (Goetze and Evans, 1979; Chappie 
and Forsyth, 1979; Bodine et. al., 1981) include a variable 
rheology with depth and presumably yield more realistic 
estimates of flexural stresses in the lithosphere. We 
should note that elastic models of lithospheric flexure have 
been applied to phenomena other than Outer Rises by Gunn 
(1943) and Walcott (1970). For a comprehensive review of 
lithospheric flexure see Turcotte (1979). 

Quite naturally these lithospheric models have been 
tested almost exclusively on bathymetric profiles. However, 
the Outer Rise also has a clear gravitational expression. 
Watts and Talwani (1974) showed that free-air gravity 
anomalies, which they referred to as the Outer Gravity High, 
have a shape and breadth which is similar to that of the 
underlying Outer Rise/slope complex. They also indicated 


1 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Chat the Outer Gravity High can be attributed largely to 
flexure of the oceanic lithosphere. Outer Gravity Highs 
have amplitudes as large as 60 mgal and wavelengths of 
several hundred kilometers. The gravity highs of course 
possess geoldal counterparts which have similar wavelengths 
and amplitudes as great as 6 meters. We shall refer to 
these as Outer Geoid Highs. The altimeter on the SEASAT 
satellite (see Tapley et. al., 1982) has made precise, 
direct observations of the marine geoid between +72° and 
-72° latitude. These altimeter data provide dense, uniform 
coverage of all Outer Rises and are well suited to the task 
of directly testing lithospheric flexure models. Direct, 
altlmetrlc observations of the Outer Geoid High also have 
the advantage that they are less affected by seamounts, 
variations In sediment thickness, etc. than are 
corresponding gravimetric and bathymetric observations. Due 
to their high accuracy and nearly global coverage, SEASAT 
data significantly augment available information concerning 
the Outer Rises. 

In this paper we will be testing the elastic plate 
model by using the SEASAT altimeter data. Bodine et al. 
(1981) have found that the model of a uniform, thin elastic 
lithosphere overlying a fluid-like as thenosphere is adequate 
for describing flexure of the oceanic lithosphere. They 
reached this conclusion after comparing elastic model 
results with those from more complex, experimentally- 
constrained lithospheric rheologies (i.e., layered 
cataclastic - elastic - plastic). However, the more complex 
rheological models do produce more realistic stress 
estimates, and they produce lower estimates, in particular, 
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of stresses associated with horlsontal compressive loads 
transverse to the trench. Based on their comparison of the 
elastic model with the complex rheological model, Boalne et 
al. (1981) were also able to postulate a relation between 
effective elastic thickness of the lithosphere and its age 
(see also Caldwell and Turcotte, 1979). We will also be 
using our results from SEASAT altimeter data to test this 
thl ckneas-age relation. 
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As indicated in the Introduction! flexure of oceanic 
lithosphere seaward of deep ocean trenches can be adequately 
modelled as an end-loaded elastic plate overlying a fluid 
asthenosphere . Using the sane elastic plate model as that 
which was first applied to the Outer Rise problem by Hanks 
(1971) one obtains the following familiar expression for 
deflection! w, of the plate 

w - w q sin (^) exp (- (1) 


where x is the horizontal coordinate normal to the trench 
axis! w is an amplitude parameter equal to 2 */^ eX p times 
v^, the outer rise height, and a is a measure of the 
flexural wavelength which can be related to rise crest 
position by a - 4x^/n (see Fig* 1)* a can also be 

related to plate thickness, T e , by 


3<l-v 2 )*<P n -P w )] l/3 4/3 

- a 

e E 


where E is Young's modulus, v is Poisson's ratio, g is the 

acceleration of gravity and p and p are mantle and 

m w 

seawater densities. The deflection solution, (1), assumes 
the lithosphere to be loaded vertically; l.e., horizontal 
loading is negligible (see Caldwell et. al., 1976). For 
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computations, we will assume; P m ■ 3.3 x 10^ kg»ra“^, - 
1.0 x 10 ^ kg‘m~^, E ■ 6.5 x 10 ^ N»m“2 t v ■ 0.25 and g * 9.8 

m # sec“2 . 


Of course this deflection of the lithosphere gives 
rise to gravity anomalies and geold undulations. 
Furthermore, surface observations of gravity (see Watts and 
Talwanl, 1974) suggest that the Outer Rise and outer trench 
slope are not isostatlcally compensated. Rather, observed 
gravity Is consistent with a model (see Figure 1) of 
uncompensated topography wherein the crust and lithosphere 
deflect as a unit. The relationship - more specifically the 
admittance - between bathymetry and geold heights derived 
from SEASAT dita Is also consistent with this model of 
uncompensated lithospheric flexure (see Appendix). 


In order to derive the geold undulations for this 
elastic plate model, two approaches have been tried. The 
first Involves numerical integration. Geold undulations 
arise due to density anomalies created by deflection of the 
seafloor and ,.he Moho. The Integral expression for these 
undulations Is 


AN(x ) 



d 

1 (p -p )lnRdz 

d-w(0 


d+t 

+ / (p -p )£nRdz 

d+t-w( 5) 


d£ + L( x ) 


(3) 
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In the above expression, d Is undeflected seafloor depth* t 

is crustal thickness* v(£) Is given by equation (l) and 

L(x) Is the contribution from mass landward of the trench 

axis (x < -x ; see discussion of mirror image assumption 
o 

below). Numerical Integration of (3) gives results that are 
quite consistent with second approach - the analytical 
approximation represented by equation (8) below. However* 
the numerical Integration approach Is not computationally 
efficient, so we have adopted a second approach. We use a 
thin mass layer approximation developed by Ockendon and 
Turcotte (1977). Deflections, w(x), of the lithosphere 
shall be assumed to manifest themselves gravitationally as 
two mass layers, the first layer resldl g at depth d and 
possessing a variable mass per unit area 

o u (x) - (P C "P W > w(x) , (4a) 

the second layer residing at Moho depth, d + t, and 
possessing a mass per unit area 

® £ ( x ) - (P m “P c ) w (x) (4b) 
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Ockendon and Turcotte (1977) represented geold undulations 
due to thin nets layers with a series expansion In powers of 
relative thickness. Retaining first and second order terns 
In the expansion (see Haxby and Turcotte, 1978), the 
following expression for geoid undulations, AN(x), results 


AN(x) - li£ | / p (l) (?) in 1 


x-C) 


d £ . . . 


( 3 ) 


- w 


p (2) (x) - P (1) (x)(d+t) 


1 ! 


In the above, p^ 2 ^(x) and p^^(x) are the equivalent source 
and dipole aonent of aass per unit area and are given by 


P 


( 1 ) 


(x) 


d+t 

/ Ap(x , z )dr 
0 


p (2> 


u> 


l+t + t' 

j j Ap(x,z)dzdz' 
0 0 


where 


(6a) 


(6b) 


A p(x , z ) - o u (x)6(z-d) + o t (x) 6(z-d-t) 


7 




and 6 la the delta function* Combining equations 1, 4, 5, 
and 6 yields 


AN(x) 





-X 


-x 


+ (p -p )w 
a w o 


1° sin 


-£-2x 


o) ( 

— ' exp v 


t+2x 


o) 


An <n 


ir w sin — exp (— ) 
o a r a 


(p c -p w ) 


t + (p. 


HI 


(7) 


In (7) the plate deflections, w(x), are assumed to be 
symmetric about the trench axis x ■ -x 0 where w(x) Is given 
by (1). This mirror Image assumption Is a far better 

approximation of reality than is the simpler approximation 
v(x) - 0 for x < -x 0 . If for simplicity one were to assume 
w(x) - 0 for x < -x 0 the second terra in (7) would be 
dropped . 


The two integrals in (7) can be evaluated exactly. 
The following expression for geold undulations results: 


AN( x ) - 


( p -p ) Gw a 
' m w 7 o 


_ 0 o 

2exp — cos — 
r a a 


sin 



a 


In 


< x +x o > 

r* 
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- exp (~)(coa i - sin 




- exp <^)<coe ~ + sin £) . Re ^ | 


(x+x 0 )(-l*»-l ) 


))••• 


♦ exp 


(2x +x)/ (2x +x) ( lx +x) 


V 


cos — - — + tin 


).i. ( El (. 


(x+x 0 )(l+i) 


I 


( 2x +x) / ( 2x +x) (2x +x) 

of o .0 

-exp — - — co* — *in — — 


) ' r, M' 


(x+x o )(l+i)' 


— w *1 n * exp<~) {(p -p )d + (p -p ) tl 

go O rN 0 * ' W ' C W 1 


(«) 


where E,(0 1* Che exponendal Integral given by 


^(O 


/ 

* 



dt 


and € is, in general, a complex variable (aee Abramovltz and 
Stegun, 1972 , p. 228) • Re( ) and Im( ) in (8) denote 
respectively the real and imaginary parts of the particular 
argument. Thla expression, (8), for geold undulations due 
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to plate flexure can be compared directly with SEASAT 
altimeter height profiles from which a long wavelength geold 
has been subtracted. This comparison can provide estimates 
of model parameters such as flexural wavelength, a (or 
equivalently plate thickness, T e , see eqn. 2) or amplitude, 
w Q . Such estimates have been obtained In the Results 
section using the method of least squares. Construction of 
the normal equations for least squares requires expressions 
for partial derivatives of N(x,w 0 , x 0 , a,t) with respect to 
parameters and t. Analytical expressions for these 
partials are obtained by differentiating (8) with respect Lo 
the parameters which are being estimated. 


10 


ORIGINAL 
°F POOR 


PAGE fs 
QUALITY 


SEASAT OBSERVATIONS 


Segments of 102 SEASAT altimeter passes over six 

trench/Outer Rise complexes have been analyzed (see Figure 
2a,b,c). For each trench the number of passes analyzed Is 
as follows: Middle America, 14; Philippine, 12; Aleutian, 

23; Mariana, 12; Izu-Bonin, 10; Kuril, 30. The sea surface 

height profiles obtained from these passes will be treated 
as geold height profiles. If the oceans were completely 

passive, the mean sea surface that SEASAT observed would, by 
definition, coincide with the geold. For this analysis, the 
ocean dynamic component of sea surface height Is, in most 
instances, the major error source. Mesoscale oceanographic 
phenomena such as ring currents, eddies, etc., constitute 

the bulk of this oceanographic "noise". Long 

wavelength (>1000km) phenomena will have no effect as such 
wavelengths are not Included In this analysis. Amplitudes 
of these mesoscale oceanographic features may exceed a meter 
in places (e.g., the Oya shl o-Kur o shi o confluence zone near 
the Kuril trench; see Cheney, 1977). In most places these 
amplitudes are substantially less than a meter. The 
precision of the actual sea surface height measurement 
after GDR corrections - is about 10 cm (Tapley et al., 
1982) . 


For comparison with the Outer Geold High model (8), 
pass segments 100 sec (700 km) in length have been chosen. 
These pass segments cover exclusively oceanic lithosphere on 
the verge of subduction; the segments terminate at the 
trench axis. These 100 sec segments have been high-pass 


11 


^oifttNAL PAGE 18 


filtered. This has been accomplished by removing a degree- 
and-order 48 background geold. The background geoid chosen 
i 8 a truncation of an updated version (Frey and Klosko, 
1983) of GEM IOC (Lerch et al., 1981) which uses Rapp's 
(1982) SEASAT - derived gravity anomalies at high degrees 
and orders. This high pass filtering, of course, removes 
unoodelled long wavelength effects including geold signals 
associated with the subducting slabs and the thickening of 
the thermal lithosphere with age. 
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RESULTS 


The 102 SEASAT pass segments of sea surface heights 
which are described in the Observations section above (see 
Figure 2a,b,c) have been fitted with the flexural model 
given by Equation (8) . The method of least squares was used 
to recover estimates of three model parameters plus a bias 
(or vertical offset) parameter for each pass. The three 

model parameters are the flexural wavelength, a, the 
amplitude, w^, and the trench axis position, x Q . Efforts 

were made to recover a fourth model parameter - crustal 
thickness or Moho depth, t; these efforts were not 
successful. This can be attributed to the fact that the 
model results (8) are rather insensitive to crustal 

thickness. 

One or two representative passes from each of the six 
trenches are shown In Figures 3a-3i ; both the SEASAT profile 
and the best-fitting model result are shown. Table I 
summarizes average values of model parameters, a and w^, 
derived from the SEASAT passes over each of the six 
regions. We will discuss the Aleutian region in detail, 

below; the other five regions will be discussed briefly. 


1. Aleutian Trench 

Twenty-three SEASAT passes across the central Aleutian 
trench (see Figure 2a) have been analyzed. The individual 
estimates of flexural wavelength, a, that were obtained 
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after least squares adjustments, are shown In Figure 4a. 
The values of o which are shown in Figure 4a have not been 
corrected for non-perpendicularity of the SEASAT groundtrack 
to the trench axis. This geometric correction is 

accomplished by multiplying times the sin$ where $ is the 
angle between the groundtrack and the trench axis. For 
Aleutian passes this geometric correction is essentially 
negligible (i.e., $~90°). Note in Figure 4a that a varies 
only slightly from pass to pass, which is consistent with 
the small "formal" errors. 

The least squares adjustment procedure was also used 
to independently recover estimates of the same three model 
parameters, a, w Q t and x Q , plus a bias parameter for the 
interpolated bathymetry profiles which coincide with the 23 
SEASAT passes. The recovered values of flexural wavelength, 
a, are plotted in Figure 4b. Although there is fair 
agreement, particularly west of the Amlia fracture zone, 
between the bathymetrically and a ltimetrically derived 
flexural wavelengths, the scatter in the 23 values obtained 
from bathymetry is significantly greater than that exhibited 
by the values derived from SEASAT. Caldwell et al. (1976) 
obtained a value of a » 67 km from a SEAMAP13 shiptrack 

bathymetry profile across the Aleutian trench. This 
shiptrack value is shown in Figure 4b for comparison with 
interpolated profiles. It is clearly consistent with 
Interpolated profiles in the vicinity. A nearby SEASAT 

profile, rev 582 (see Figure 3a), yields an estimate 
of a ■ 91.8 km which is significantly larger than Caldwell 
et al.'s (1976) estimate from bathymetry. In general 
estimates derived from bathymetry are significantly less 
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consistent from pass to pass than are those derived from 
SEASAT. Furthermore, the estimates of a which we obtained 
from bathymetric fits have a tendency - to be less than 
those derived from altimetry. This tendency Is exhibited at 
other trenches and is discussed below. Seamounts, sediments 
and other localized topographic variations act as noise 
obscuring the Outer Rise; these same noise sources are 
significantly attenuated in the geoid signal, i.e., the 
Outer Geoid High. In some places, as discussed below, an 
Outer Geoid High is evident In the SEASAT data where an 
Outer Rise can not be found in the bathymetry. 



As 

an 

indication of the good 

agreement 

between 

the 

model 

and 

the 

SEASAT data, It should 

be noted 

that the 

rms 

diffe 

r enc-j 

bei 

tween model and data ra 

nges from 

22 to 40 

cm 


for the Aleutian passes. For rev 1185 (see Figure 3b) - 
the westernmost of the 23 - the ms difference Is 22 cm 
which is not much larger than the measurement precision of 
SEASAT. 

The Amlia fracture zone Intersects, obliquely, the 
Aleutian trench between 186° and 187° E. This fracture zone 
- and possibly its down-slab extension - could be acting as 
a line of weakness in the lithosphere. It appears, in any 
case, to be affecting the Outer Rise and the Outer Geoid 
High. Therefore, the three SEASAT passes, revs 654, 166 and 
1156 (dashed in Figure 2a), which intersect the trench 
between 186° and 187° E are excluded from the ensemble of 20 
passes used to construct mean values of model 
parameters, o and w Q > In Table I. 


IS 
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Thirty ascending passes of SEASAT data across the 
Kuril trench were analyzed. However, the eight 

southwes ternnos t. passes of these 30 were culled before 
regional mean values of a and w Q were computed for Table 
I. These eight parses cover (see dashed groundtracks in 
Figure 2a) Outer Rise which is influenced by subduction at 
the Japan trench as well as the Kuril trench. The Japan 

trench axis is oriented at roughly a 45° angle to the 
groundtracks of these passes - whereas the Kuril trench is 
nearly perpendicular to them. The high values (range 118 to 
136 km) of flexural wavelength, a^, (uncorrected for 
groundtrack geometry) are probably due to the influence of 
the Japan trench. It is also conceivable that high 

mesoscale variability in this region (see Observations 
section above, Cheney, 1977) might also be affecting these 
passes . 


Recovered values of flexural wavelength, a, for the 
remaining 22 SEASAT passes across the Kuril trench range 
from 88 to 106 km with a mean of a ■ 96 km. Using 10 

shiptracks of bathymetric data, McAdoo et al. (1978) 
recovered values which range from 67 to 102 km with a mean 
of 89 km. As in the case of the Aleutians, the 

altimetrically derived estimates of a are more consistent 
from pass to pass than are the bathymetrically derived 
estimates. Furthermore, the maximum values recovered for 
both the Kuril and the Aleutian regions from SEASAT data 
nearly equal those from bathymetry. Once again, it appears 
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then that bathymetry tends to produce low estimates of 
flexural wavelength a. We will discuss this further, 
below. Figure 3c shows a representative pass, rev 469, 
across the Kuril trench. For this pass a is estimated to be 
102 km and the rms difference between model and SEASAT 
profile 1 8 52 cm; this rms difference ranges between 20 and 
55 cm for the 22 passes. 


3. Middle America Trench 

Passes on or northwest of the Tehuantepec ridge were 
also analyzed but discarded Inasmuch as these passes 
traversed the crest of the East Pacific Rise and produced 
results that were not consistent from one pass to the 
next. On the other hand, all descending SEASAT passes 
between the Tehuantepec and Cocos ridges (see Figure 2b), 
excluding redundant coincident passes, were analyzed and 
retained. These passes were fourteen in number. 

These fourteen passes produce quite consistent 
results. The mean flexural wavelength, a, (before slight 
correction for groundtrack geometry) is 73.8 km for this 
sample of fourteen; the standard deviation Is 4.4 km. Rev 
558 is shown in Figure 3d as a representative of sea surrace 
and model geold heights. The rms difference between model 
and observation for 558 Is 28 cm; the rms differences for 
the rest of the passes range between 16 and 47 km. It 
should also be noted that we were not able to estimate 
flexural wavelengths for the interpolated bathymetric 
profiles along many of these 14 groundtracks as they did not 
even display an Outer Rise. 
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4. Philippine Trench 



Eight SEASAT passes which traverse the Outer Rise 
seaward of the Philippine trench were analyzed. The 

resulting estimates of flexural wavelength were less 
consistent from pass to pass than similar estimates for 
other trenches. The mean value of flexural wavelength 

(uncorrected for groundtrack geometry) Is 120.0 km and the 

standard deviation is 14.5 km. These passes are oriented at 
approximately a 34° angle with respect to the trench axis 
(see Figure 2c); therefore, the geometric correction is 
large. Two SEASAT profiles, rev 578 and 736, are shown in 
Figures 3e and 3f. The agreement between model and 

observation is particularly good for rev 578; the rms 
difference is 13 cm. It is interesting to note that 

coincident profiles of interpolated bathymetry yield poor 
results; some passes exhibited no Outer Rise and the 
recovered flexural wavelength estimates, (uncorrected 

for geometry), range from 50 to 127 km. 


5. Izu-Bonin Trench 

The presence of the Ogasawara Plateau and Ramapo Bank 
apparently interfere with development of a typical Outer 
Rise and Outer Geold High adjacent to the southern Izu-Bonin 
trench. Bodine and Watts (1979) found that it is difficult 
to resolve an Outer Rise in shiptrack bathymetric profiles 
across this area. Only four SEASAT passes could reasonably 
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be analyzed and theae were cloaely spaced. One of these, 
rev 664, is shown In Figure 3g. Note that in part because 
of the acute angle (32°; see Figure 2c) between groundtrack 
and trench axle, the Outer Geold High la quite long. For 
the four passes, the mean value of uncorrected flexural 
wavelength, a u , Is 237.7 km. This uncharacteristically 
large value of "a u (due to geometry plus an apparently thick 
lithosphere) necessitated using longer (128 sec) pass 
segments for this region. The corrected mean flexural 
wavelength Is a • 125.9 km. 

Bodlne and Watts (1979) found that It la difficult to 
estimate flexural wavelengths from shiptrack profiles of 
bathymetry across the Izu-Bonin trench. They handled this 
difficulty by assuming, a priori, a flexural 
wavelength, a (or x^), which is consistent with a relation 
between age and effective lithospheric thickness. This 
a priori flexural wavelength is 90 km — considerably less 
than our estimate from altimetry (126 km). 


6. Northern Mariana Trench 

Nine descending passes of SEASAT data across the 
northern Mariana trench were analyzed (see Figure 2c). 
Consistency of flexural wavelength estimates from pass to 
pass is poor. Because of the curvature of the trench, 
geometric corrections had to be applied Individually to each 
pass. Due to a very acute angle between groundtrack and 
trench axis, three passes, revs 721, 477 and 233 (see dashed 
groundtracks in Figure 2c), were culled from an original 
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group of 12. A mean flexural wavelength (corrected for 
geometry) of 97.9 km (o - 22.9 km) waa eatimated ualng nine 
paasea. Interestingly, Bodine and Watts (1979) found that 
an Outer Rise was quite difficult to detect in shlptrack 
profiles of bathymetry across the northern Mariana. Two 
passes, revs 879 and 836, are shown in Figures 3h and 31 
respectively. Rev 879 agrees quite well with the 

model (<* u - 128 km; rms of fit, 35 cm); rev 836 deviates 

substantially from the model ( • 107 km; rma of fit, 94 
cm) due to the effect of two large seamounts in the western 
Marcus ridge. 
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INTERPRETATION 


Bodlne et *1. (1981) applied a nore realistic, complex 
rheological model of oceanic lithosphere to the phenomena of 
the Outer Rise. They resorted to a three-layer rheology - 
the shallowest layer being cataclastlc, the next elastic and 
the deepest plastic. Rock mechanics experiments by Goetze 
and Evans (1979) on olivine were used to constrain 
characteristics of the plastic layer. Bodlne et al. (1981) 
found that the Outer Rise shape predicted by this complex 
rheological model was not significantly different than that 
predicted by the simple elastic plate model with similar 
plate thickness. In fact, they were able to use the complex 
model to predict how the effective elastic thickness of 
oceanic lithosphere should depend on temperature and hence 
age. Using the model of Parsons and Sclater (1977) to 
describe the thermal evolution of the lithosphere, Bodlne et 
al. ( 1981 ) postulated that the effective elastic thickness 
of the lithosphere, T e , Is given by 

T fi - C . age 1/2 (9) 


where C ■ 4.6 ± 0.6 km»my -1 / 2 . This thickness law, (9), was 
derived by using the results of Goetze and Evans (1979) from 
experiments on dry olivine. A wet dunlte rheology led to a 
slightly lower value of C. The thickness law based on a dry 
olivine rheology Is plotted in Figure 5. 
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For a simpler Interpretation of the thickness-age 
lav, one can think of the effective elastic thickness as 
simply the depth to a particular lsothera. Caldwell and 
Turcotte (1979) find the 700°C lsothera aost suitable. 

Our eatlaatea of effective elastic thickness derived 
fron SEASAT altimeter data are also plotted In Figure 5. 
These six estimates correspond to the six regions which we 
studied — the Kuril, Aleutian, Middle Aaerlca, Philippine, 
Mariana and Izu-Bonln. The estimates are also listed In 
Table I. They are derived froa the estlaates of mean 
regional flexural wavelength, a, by using equation (2). 

The agreement between our results for these six areas 
and the thickness age relation, (9), is quite good. The 
agreement Is In contrast to the rather rough fit obtained 
using estimates derived from bathymetric data. This rough 
fit as depicted in Figure 10 of Bodlne et al. (1981) is 
mostly a result of the low estimates of a and T # obtained 
from bathymetric profiles across the Kuril, Izu-Bonln and 
Mariana trenches. Recall that we Indicated in the Results 
Section that bathymetry, in contrast to altlaetry has an 
apparent tendency to yield low estimates of T ft o* a. Bodlne 
et al. (1981) explained these low estlaates of T e froa 
bathymetry by postulating that horizontal compression Is 
acting across these trenches (Kuril, Izu-Bonln and Mariana); 
they suggested that this compression (1 kbar) reduces the 
effective thickness. They use the complex rheological 
model - as opposed to the simple elastic model - to describe 
this reduction In effective thickness. The complex 

rheological model Is also vety useful for producing 
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realistic estimates of stress and explaining certain steep 
outer trench slopes. If our thickness estimates, T # , from 
altimetry as depicted in Figure S are valid, then it is not 
necessary to invoke horlcontal compressive loads of this 
magnitude (1 kbar) across the three trenches* 
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SEASAT altimeter data over Outer Rises apparently 
provide better estimates of effective lithospheric thickness 
than do bathymetric data. We base this observation on the 
fact that altlmetrlc results in a given region are more 
consistent from profile to profile. We have also suggested 
that bathymetric profiles across Outer Rises have a tendency 
to occasionally produce values of flexural wavelength - or 
plate thickness which are significantly less than 
corresponding estimates from SEASAT data. 

We can speculate about reasons for this apparent 
superiority of altimeter data as an estimator. Altimeter 
data enjoys the advantage of the geoldal representation in 
which unrodeled, high frequency effects such as seamounts 
are attenuated. At present, however, we do not fully 
understand this discrepancy between altlmetrlc and 
bathymetric estimates of plate thickness. 

Regarding bathymetric profiles across Outer Rises, 
Jones et al. (1978) showed that the confidence limits on 
estimates of tlexural (or plate thickness) are, in some 
cases, quite broad — leading to uncertain estimates. Agnew 
and Sharman (personal communication, 1983) found that this 
uncertainty is sufficient in certain cases to permit only 
estimation of a lower bound on plate thickness. 
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CONCLUSIONS 


SEASAT altimeter data provide significant new 
Information on the deformation of the lithosphere seaward of 
deep ocean trenches. These data, when compared with a 
mechanical model, produces estimates of the effective 
elastic thickness of the lithosphere which are more 
consistent from profile to profile in a given region than 
those derived from bathymetry. For the six trenches which 
we studied, the effective elastic thickness derived from 
altimetry increases with the square root of lithospheric age 
as predicted by Bodine et al. (1981). Significant 
horizontal compression across these trenches is not 
indicated. This is consistent with Forsyth's (1982) studies 
of Outer Rise earthquakes from which he concluded that 
regional compressive stresses must be a small fraction of 
the bending stresses. 
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APPENDIX 

1. A TEST OF THE MODEL OF UNCOMPENSATED TOPOGRAPHY 

Watts and Talvanl (1974) proposed a model to explain 
the relationship between gravity anomalies and topography 
seaward of deep-ocean trenches. This model which they 
referred to as the crustal upwarplng model (their #3) is 
essentially the same model as that which we are using in 
this study (see Figure 1); i.e., the crust and the 
lithosphere flex as a unit. We wish to test whether this 
model also can describe the observed relationship between 
geold undulations observed from SEASAT and topography. We 
will represent this relationship in wave number (k) domain 
by using the admittance function, Z(k); this representation 
has been popularized for studies of gravity by McKenzie and 
Bowln (1976) and others. Geoid undulations and bathymetry 
are characterized, respectively, as output and input signals 
in a linear system. The spatial Fourier transform of 
geold, A(N), is given by 


F ( AN( x ) ) - Z(k) . F( B(x ) ) 


(A-l) 


where F(B(x)) is the Fourier transform of bathymetry. The 
theoretical admittance function for our crustal upwarplng 
model (Figure 1) is 


.L ' V. ' U ,\| /i' JJLME77 
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Z(k) - e- kd l(e c -P B ) + (p.-p,,) e' kt ] (A-2) 


which can be obtained by dividing McKenzie and Bowin's 
(1976) expression (number 14 for admittance between gravity 
anomalies and bathymetry) by g k The theoretical 

admittance is plotted as a continuous curve in Figure 6. \ 

crustal thickness, t, of 10 km is assumed; all other 
parameters have the values identical to those which were 
used In the Model section above. 

Figure 6 also shows amplitude estimates of admittance 
between SEASAT observations of sea surface height and ocean 
depths Interpolated from NORDA digitized bathymetry. These 
estimated admittances are obtained from 23 SEASAT passes 
over the Aleutian trench (see Figure 2a). For each pass a 
sequence of the 128 (2^) data points just seaward of the 
trench axis Is formed. A mean was removed from each 
sequence and a degr ee-and-order 12 background geoid was 
subtracted from each height profile. The 23 sequences were 
then combined Into a single sequence of SEASAT data; a 
corresponding bathymetric sequence was also formed. After 
the fashion of McKenzie and Bowln (1976) admittance 
estimates were obtained from the ratio of cross-spectral 
estimates to the bathymetric power spectral estimates. To 
obtain estimates of the cross spectrum, Welch's method of 
averaged modified periodogram was used In conjunction with 
fast Fourier transform algorithms. 
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The estimated admittances in Figure 6 agree 
reasonably well with theoretical admittances thereby 
supporting the popular assumption that topography of the 
Outer Rise Is uncompensated. The crustal upwarplng model Is 
consistent with the admittances estimated from SEASAT 
data. Admittances estimated from altimeter data over mid- 
plate (compensated) topography (Watts, 1979) are much less 
at wavelengths In the range 300 to 1000 km. 
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TABLE I 


Trench 


Flexural 

Wavelength 

ot(km) 

Amplitude 

V“> 

Elastic 

Thickness 

T e (km) 

Age 

(10 6 Yr) 

Mariana (9) 


107.9 

465.+ 75. 

51.3+0.3 

140-160 

Kuril (22) 


92.7 

513. ±139. 

41.9+ 5.1 

100-110 

Aleutian (20) 

90.8 

461. ±61. 

35.4+ 1.4 

52- 65 

Philippines 

(8) 

66.3 

307 .±74. 

26.8+ 4.3 

50- 55 

Middle Amer 

(14) 

71.7 

322 .±46. 

29.7+ 2.4 

22- 30 

Bonin (4) 


126.0 

354 .±20. 

63.0+ 8.6 

135-160 
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Figure 2b. Qroundtrectei o I tCAIAT pwm gmlywd owr the Middle Araertoe trench, lee text 
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MARIANA TRENCH 
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Figure 2c. Groundtrackt of 8EA8AT pane* analyzed over the Philippine, Izu Bonin end Marlene tranche*. 

See caption for Figure 2a. 
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Figure 3a- 3b. Selected SEASAT passes over trenches indicated. Residual sea surface heights 
(with degree-and-order 48 background geoid removed) shown as solid curves. Best-'itting 
theoretical geoid profiles from elastic plate model are shown as dotted curves. 
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Figure 3c-3d. Selected SEASAT pet set over trencher indicated. Residual tea surface heights 
(with degree-end-order 48 background geo id removed) shown at solid curves. Best-fitting 
theoretical geoid profiles from elastic plate model are shown as dotted curves. 
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Figure 3e-3f. Selected SEASAT passes over trenches indicated. Residual see surface heights 
(with degree-and-order 48 background geoid removed) shown as solid curves. Best-fitting 
theoretical geoid profiles from elastic plate model are shown as dotted curves. 
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Figure 4a. Value* of flexural wavelength, o^, estimated for each of 23 
SEASAT passe* over the Aleutian trench plotted vs longitude of 
intersection between trench axis and ground track. 



178 ° 180 ° 182 ° 184 ° 186 ° 188 ° 190 ° 192 ° 194 ° 

EAST LONGITUDE OF TRENCH INTERSECTION 
DEGREES 


Figure 4b. Values of flexural wavelength, a u , estimated for 23 interpolated 
bathymetric profiles coinciding with 23 SEASAT passes used in Figure 4a. 
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